Magnetic degeneracy and hidden metallicity of the spin density wave state in 

ferropnictides 
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We analyze spin density wave (SDW) order in iron-based superconductors and electronic structure 
in the SDW phase. We consider an itinerant model for Fe— pnictides with two hole bands centered 
at (0,0) and two electron bands centered at (0, tt) and (it,0) in the unfolded BZ. A SDW order 
in such a model is generally a combination of two components with momenta (0,7r) and (tt, 0), 
both yield (tt, tt) order in the folded zone. Neutron experiments, however, indicate that only one 
component is present. We show that (0,7r) or (tt, 0) order is selected if we assume that only one 
hole band is involved in the SDW mixing with electron bands. A SDW order in such 3-band 
model is highly degenerate for a perfect nesting and hole-electron interaction only, but we show that 
ellipticity of electron pockets and interactions between electron bands break the degeneracy and 
favor the desired (0, tt) or (tt, 0) order. We further show that stripe-ordered system remains a metal 
for arbitrary coupling. We analyze electronic structure for parameters relevant to the pnictides and 
argue that the resulting electronic structure is in good agreement with ARPES experiments. We 
discuss the differences between our model and Ji — J2 model of localized spins. 

PACS numbers: 74.70.Xa, 75.30.Fv, 75.25.-j 



I. INTRODUCTION 

The discovery of superconductivity in the oxypnictide 
LaFeAsO (Ref. [0) created a new class of Fc-based high- 
T c superconductors - ferropnictides (FPs) The phase 
diagram of FPs is similar to that of high-T c cuprates 
and contains an antiferromagnetic phase at small dop- 
ings and superconducting phase at larger dopings. There 
are two important distinctions, however. First, parent 
compounds of FPs are antiferromagnetic metals, and sec- 
ond, the pairing symmetry in FPs is, most likely, an ex- 
tended s— wave, with or without nodes. — ^ Electronic 
structure of parent FPs in the normal state has been mea- 
sured by angle-resolved photoemission (ARPES)^ and by 
magneto-oscillations^. It consists of two quasi-2D near- 
circular hole pockets of non-equal size, centered around 
r point ((0,0)), and two quasi-2D elliptic electron pock- 
ets centered around (0, ±7r) and (±7r, 0) points in the 
unfolded Brillouinc zone (BZ) which includes only Fe 
atoms. For tetragonal symmetry, the two electron pock- 
ets transform into each other under a rotation by 90°. In 
the folded BZ, which is used for experimental measure- 
ments because of two non-equivalent As positions with 
respect to an Fe plane, both electron pockets are cen- 
tered around (7r, 7t). The dispersions near electron pock- 
ets and near hole pockets are reasonably close to each 
other apart from the sign change, i.e., there is a substan- 
tial degree of nesting between hole and electron bands. 
There also exists the 5-th hole band near (7r,7r), but it 
is more three-dimensional and does not seem to play a 
significant role, at least for magnetism. 

In this paper we analyze spin and electronic struc- 
tures of a magnetically ordered state in parent FPs below 



Tsdw ~ 150Ki£i We will only focus on FeAs mate- 
rials. Neutron scattering measurements on parent FeAs 
pnictides have found that the ordered momentum in the 
unfolded BZ is either (0,7r) or (tt,0), i.e., magnetic or- 
der consists of ferromagnetic stripes along one crystallo- 
graphic direction in an Fe plane, and antiferromagnetic 
stripes along the other direction. 

Such magnetic order emerges in the J\ — J2 model 
of localized spins with interactions between nearest and 
next-nearest neighbors (J\ and J2, respectively), for 
J 2 > 0.5Ji4>2. A localized spin model, however, is best 
suitable for an insulator and is generally not applica- 
ble to a metal unless the system is close to a metal- 
insulator transition.— An alternative scenario, which we 
explore here, assumes that parent FeAs FPs are "good 
metals" made of itinerant electrons, and antiferromag- 
netic order is of spin-density-wave (SDW) type. Given 
the electronic structure of FPs, it is natural to as- 
sume that SDW order emerges, at least partly, due to 
near-nesting between the dispersions of holes and elec- 
tron o 10 ^ 1 ! 12 : 13 ! 14 : 15 : 16 : 17 . Such nesting is known to give 
rise to magnetism in Cr (Ref. OH). The itinerant scenario 
for FPs is largely supported by a reasonable agreement 
between ab initio electronic structure calculations 1 ^ and 
magnetooscillation and ARPES experiment a 20 i 21 ' 22 i 23 ' 24 , 
although electronic structure calculations also indicate 
that magnetism partly comes from other fcrmionic bands 
which do not cross the Fermi level^. 

We discuss below several puzzling electronic and 
magnetic features of the SDW state which have not 
been yet understood within the itinerant scenario. They 
include: 
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• A higher conductivity in the SDW state - the sys- 
tem becomes more metallic below Tsdw- Naively, 
one would expect a smaller conductivity due to at 
least partial gapping of the Fermi surfaces (FS). 
A reduction of the quasiparticle damping, induced 
by such gapping, may slow down the decrease of 
conductivity, but the observed increase is highly 
unlikely for a conventional SDW scenario. 

• A significant reconstruction of the electronic dis- 
persion across the SDW transition and compli- 
cated FS topology below Tsdw, with more FS 
crossings than in the normal state (Refsi 28 i 29 ' 30 i 31 ) 
Some ARPES measurements sho w 28 ' 29 i 30 that in 
the SDW state visible dispersion remains hole-like 
around the T— point despite apparent mixing of 
electron and hole bands, while close to M = (tt, tt) 
both electron and hole bands are present. Other 
measurements 3 -! indicate that hole and electron 
pockets are present both near M and near T points. 

• The particular (0, tt) ((tt, 0)) ordering of Fe spins 
in the unfolded BZ. It is straightforward to ob- 
tain (tt, tt) ordering in the folded BZ simply be- 
cause hole and electron bands are shifted by (tt, tt) 
in the folded zone. However, SDW order in the 
folded zone involves two separate magnetic sublat- 
tices, and unfolding reveals two SDW order pa- 
rameters (OP) Ai and A2 with momenta = 
(0, tt) and Q2 = (tt, 0), respectively (see Fig. [1]) 
A generic spin configuration then has the form 
S(R) = AV QlR + A 2 e^ 2R (the sublattice OPs 
are Ai + A2 and Ai — A2). For such configura- 
tion, next-nearest neighbors are antiferromagneti- 
cally oriented, but the length and the direction of 
S between nearest neighbors varies. Only if either 
Ai or A2 vanish, OPs of the two sublattices align 
parallel or antiparallcl to each other, and the spin 
configuration becomes the same as in the experi- 
ments. The issue then is what interaction causes 
Ai or A2 to vanish. 

Another unsolved issues which we do not address here 
are the absence of magnetism in the nominally undoped 
LaFePO and LiFeAs compounds and the type of a mag- 
netic order in Fei+ x Te/Se systems. The absence of mag- 
netism in LaFePO and LiFeAs could be due to the fact 
that these materials are less quasi-two-dimensional, with 
less nesting^ For Fei +a; Te/Se, the experimental data are 
still controversial^ 

The full analysis of possible SDW ordcrings in FeAs 
systems within a generic 4-band model is quite messy, 
and one should have a good staring point to be able to un- 
derstand the physics. One way to analyze the problem is 
to solve first the model of 4 equivalent 2D circular bands 
(two hole bands and two electron bands) with isotropic 
interband and intraband interactions, and then introduce 
anisotropy between the two hole bands, ellipticity of the 
two electron bands, and anisotropy of the interactions. 



This approach has been put forward by Cvetkovic and 
Tesanovio 3 - 2 -. They solved exactly the isotropic model 
and found that SDW state is an insulator - all 4 Fermi 
surfaces are fully gapped, and that there exists a degener- 
ate manifold of SDW orders. The degeneracy is the same 
as in J2 model of localized spins: the SDW order involves 
two antifcrromagnctic sublattices with equal magnitudes 
of the OPs (i.e., Ai ■ A2 = 0), but the angle between 
the two sublattices (j) = cos" 1 (A 2 — A^)/(A 2 + A|) can 
be arbitrary. The (0, tt) and (tt, 0) states belong to this 
manifold, but are not yet selected in the isotropic model. 
The issue not yet addressed within this approach is what 
happens with a degenerate manifold once one moves away 
from the isotropic limit. 

Our starting point is different. We use the experi- 
mental fact that the two hole FSs are of quite differ- 
ent sizes 2 ^ and assume that one hole band interacts with 
electron bands much stronger than the other. We then 
consider, as a first approximation, a 3-band model of 
one hole band centered around T-point and two electron 
bands centered around (0,7r) and (tt, 0). We solve this 
model in the mean-field approximation and show that 
the SDW order is still degenerate for isotropic interac- 
tions and circular bands. A degenerate manifold consists 
of two antifcrromagnctic sublattices with generally non- 
equal magnitudes of the sublattice OPs. The manifold 
includes (0, tt) and (tt, 0) states among many others. 

In terms of Ai and A2, the degeneracy implies that 
the ground state energy depends only on the combina- 
tion (Ai) 2 + (A2) 2 . The excitation spectrum of such de- 
generate state contains five Goldstone modes. We then 
show that the ellipticity of the electron bands and the 
anisotropy of the interactions breaks the degeneracy and 
adds to the energy the term /3i2|Ai| 2 |A2| 2 with posi- 
tive (3i2- The full energy is then minimized when either 
Ai = 0, or A2 = 0, ellipticity and the anisotropy of the 
interactions both select (0, tt) or (tt, 0) stripe order already 
at a mean-field level. This is one of the main results of 
this paper. A stripe SDW order mixes the hole band 
and one of the two electron bands, but leaves another 
electron band intact. As a result, the system remains a 
metal even at strong SDW coupling. We verified that the 
selection of the stripe order by the interactions between 
electron states holds only if the interactions are in the 
charge channel. The same interactions, but in the spin 
channel, select a different state with Ai = iAa, in which 
SDW order resides in only one of the two sublattices, see 
Fig. [3d). 

We did not consider corrections to the mean-field the- 
ory (i.e., quantum fluctuations). They are relatively 
small, at least at weak coupling, although, very likely, 
they also break 0(6) symmetry. We also do not consider 
the coupling to the phonons, which is another poten- 
tial source of symmetry breaking and also pre-emptive 
tetragonal to orthorhombic transition^ 3 -. 

We next consider the role of the second hole band. This 
band interacts with the electron band left out of the pri- 
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mary SDW mixing. These two bands are less nested, and 
the interaction must exceed a threshold for an additional 
SDW order to appear. If this is the case and an addi- 
tional SDW order is strong enough, it gaps the remaining 
two FSs, and the system becomes an insulator. This sec- 
ond (weaker) SDW OP has the momentum (0, tt) if the 
original SDW order was with (tt,0), and vise versa. It 
then introduces A 2 (or Ai) which was set to zero by the 
initial selection of (0, tt) or (tt, 0) order. We show that 
this second order parameter is directed orthogonal to the 
primary one, i.e., Ai • A2 = 0. The magnitude of the 
OP at each Fe cite is then the same, but the OP is not 
ferromagnetic along x or y direction, i.e., the resulting 
SDW order is different from (0, tt) oi(tt, 0). The only way 
to preserve (0, tt) or (tt, 0) order in the 4-band model is to 
assume that the interaction between the second hole FS 
and the electron FS left out in primary SDW selection is 
below the threshold. Then A 2 (or Ai) remains zero, and 
(0, tt) or (tt, 0) order survives. A simple but important 
consequence of this observation is that the 4-band model 
with a stripe order must remain a metal. Specifically, 
one hole and one electron band are not involved in the 
SDW mixing and should be observed in ARPES exper- 
iments exactly at the same positions as in the normal 
state: the hole band near (0, 0) and the electron band 
near (tt,tt). This somewhat unexpected result is another 
main conclusion of this paper. 

Applying the results to the pnictides, where the inter- 
actions are believed to be moderate, we show that the 
hole band with a larger FS interacts more strongly with 
elliptic electronic bands than the hole band with a smaller 
FS. Hence, in our theory, the smaller hole FS stays in- 
tact below a SDW transition, and the reconstruction of 
the fermionic structure below T$dw involves the larger 
hole FS and one of two electron FSs. Since the hole 
FS is a circle and the electron FS are ellipses, the two 
cross at four points when shifted by, e.g., (0,tt), and for 
moderately-strong interactions SDW gaps open up only 
around crossing points^. In the folded BZ this gives rise 
to two FS crossings near either (0,0) or (tt, tt). Adding 
the FSs which are not involved in the SDW mixing, we 
find that there generally must be 3 FS crossings near 
(0, 0) or (tt, tt) - more than in the normal state. 

Further, our conclusion that (0, tt) SDW order in the 
4-band model preserves metallic behavior even at strong 
coupling implies that there is no continuous evolution be- 
tween our 4-band model and J\ — Ji model of localized 
spins, despite that quantum fluctuations in J\ — J — 2 
model select the same (0, it) order. At a first glance, 
this is surprising as it is well-known that in the one-band 
Hubbard model, there is an evolution from itinerant to 
localized behavior— by which we mean that, upon in- 
creasing U, a one-band system evolves from an antifer- 
romagnetic metal to an antiferromagnetic insulator with 
the same with (it, tt) magnetic order. To understand why 
our 4-band model is special, we considered a half-filled 
t — t' — U Hubbard model on a square lattice with the 
hopping t between nearest neighbors and hopping t' be- 



tween next-nearest neighbors along the diagonals. In the 
large U limit, the model reduces to J\ — J2 model with 
Ji ~ t 2 /U and J 2 ~ (t') 2 /U. As small U, the model de- 
scribes two hole bands centered at (0,0) and (tt,tt), and 
two electron bands centered at (tt, 0) and (0, tt). At t = 0, 
all four bands are identical (up to an overall sign), and 
all have an isotropic quadratic dispersion near the top of 
the hole bands and the bottom of the electron bands. For 
t 0, electron bands become elliptical near the bottom, 
and the masses of the two hole bands become different 
(i.e., one hole FS becomes larger and the other becomes 
smaller). This very much resembles the geometry of our 
4-band model, except that in t — t' — U model the second 
hole band is located at (tt,tt) (see Fig. 5). We analyzed 
t — t' — U model in the same way as the 4-band model 
and found that in t — t' — U model there is an evolution 
from a metallic to an insulating behavior as U increases. 
Namely, the state with (0, tt) or (tt, 0) stripe order is a 
metal at small U and an insulator at large U. 

The difference between our case and t — t' — U model 
can be easily understood. If we apply the same logics as 
in our case, i.e., consider the 3-band model first, we find 
(0, tt) SDW order which mixes one hole and one elec- 
tron FSs. The remaining electron band can again mix 
with the other hole band and create another SDW order, 
what eventually gaps all 4 FSs. However, in t — t' — U 
model, the second hole FS is located at (tt, tt) rather than 
at the r point, and the second SDW OP has the same 
momentum (0, tt) as the primary SDW OP. As a result, 
the two OPs just add up, and (0, tt) order survives. From 
this perspective, the reason why (0, tt) ordered state in 
the pnictides is a metal is the FS geometry: the fact that 
both hole FSs are centered at the same T point. Were 
they centered at (0, 0) and (tt, tt) the system would be 
a (0,7r) metal at weak coupling and (0,tt) insulator, de- 
scribed by Ji — J2 model, at strong coupling. 

In the rest of the paper we discuss these issues in full 
detail. The structure of the presentation is the following. 
In Sec. Ill AI we introduce three-band model and discuss 
SDW ordering in general terms. In Sec. Ill B I we consider 
the case of circular electron pockets and only the interac- 
tions between hole and electron pockets. We show that, 
at this level, there exists a degenerate manifold of SDW 
ordered states. In Sec. Ill CI we show that ellipticity of 
electron pockets and interactions between them remove 
the degeneracy and select a stripe SDW order already at 
the mean- field level. In Sec. Ill C 31 we consider spin com- 
ponent of the interaction between electron pockets and 
show that it selects a different state in which SDW order 
appears only at a half of Fe sites. 

In Sec. IIIII we add the second hole band and discuss a 
potential SDW mixing between this band and the elec- 
tron band left out of SDW mixing in the 3 band model. 
We argue that (0,tt) or (tt, 0) order survives only if this 
this additional SDW mixing does not happen, i.e., the in- 
teraction between these two bands is below the threshold. 
In Sec. IIII Al wc contrast this behavior with the one in a 
t — t' — U Hubbard model with nearest and next-nearest 
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FIG. 1: Fermi surface of ferropnictides in the unfolded BZ 
consisting of two hole pockets centered around the T— point 
and two electron pockets centered around (n, 0) and (0, n) 
points, respectively. The wave vectors Qi and Q2 are two 
degenerate nesting wave vectors. 



neighbors. We show that in this model the second hole 
FS is centered at (tt,tt) rather than at (0,0), and the 
secondary SDW OP has the same momentum as the pri- 
mary OP, i.e., the stripe order is preserved. We show 
that (0, 7r) or (ir, 0) state in t — t 1 — U model can be either 
a metal or an insulator, described at strong coupling by 
J\ — J2 model. In Sec. IIVI we consider electronic struc- 
ture of the SDW phase with (0, it) order and compare our 
theoretical results with the experiments. We argue that 
the agreement with ARPES is quite good. We present 
our conclusions in Sec. [V] 

Our results arc complimentary to earlier arguments as 
to why SDW phase remains a metal. Ran et al£§- con- 
sidered an orbital model and argued that the SDW gap 
must at least have nodes because SDW coupling vanishes 
along particular directions in k— space. Cvetkovic and 
Tesanovioi^ and Vorontsov et alM- argued that a metal 
survives even for the case of circular FSs, if SDW order 
is incommensurate. As we said, we argue that the sys- 
tem described by 4-band model remains a metal even if 
the order is commensurate (0, ir) or (ir, 0), and the inter- 
actions are angle-independent. Our analysis also serves 
as a justification to the work a 38 i 39 which discussed the 
consequences of the (n, 0) SDW order in itinerant mod- 
els without analyzing why this order is selected. Some of 
our results are also consistent with the analysis of possi- 
ble magnetically ordered states in the two-orbital model 
by Lorenzana et al^. We also found, for a particular 
model, a state with magnetic order residing only on one 
of the two sublattices. Such state was identified as a pos- 
sible candidate in a generic Landau-theory analysis for a 
two-sublattice order paramter4S 



II. ITINERANT SDW ORDER IN THE 3 BAND 
MODEL 

A. General consideration 

Consider a model of interacting fermions with a circu- 
lar hole FS centered around T-point (a-band) and two 
elliptical electron Fermi surface pockets centered around 
(±7r, 0) and (0, ±7r) points in the unfolded BZ (/3-bands) 
(see FigHJ): 
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Here, Ep 1 = 

^2 2 ^22 

tt~^~ + -n — - — u are the dispersions of hole and electron 
bands. The momenta of a— fermions are counted from 
(0, 0), the momenta of /3i— and fii— fermions are counted 
from (0, 7r) and (71", 0), respectively 

The interacting part of the Hamiltonian contains 
density-density interactions involving a and (3 fermions 
(interactions with small momentum transfer), (7r, 0), 
(0, 7r), and (n, n) scattering processes, and umklapp pair 
hopping termsi^. Consider first only the interactions be- 
tween hole and electron bands, which give rise to a SDW 
order. These are density-density interaction between a 
and P fermions, and the pair hopping terrn^: 
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We neglect potential angular dependencies oiU\ and L/3 
along the FSs. 

Because U% and t/3 involve (3\ and (3% fermions, we have 
to introduce two SDW OPs Ai oc X^pi^ipa/^ipT^) with 
momentum Qi = (0, if) and A 2 oc J2p( a ipsP2p~<&5~f) 
with momentum Q2 = (n, 0). Without loss of gener- 
ality we can set Ai along z-axis and A2 in the xz-planc. 
In explicit form, we introduce 



(3) 



AJ = -tWX>I pT /3ipT> 
p 

A z 2 {x) = -^swX^ipt^PTU))- 
p 

where Usdw = U\ +U3. 

B. Degeneracy of the SDW order 

Assume first that electron pockets are circular, i.e. 



and 



In this situation, the 
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Hamiltonian, Eq. ([2]). can be easily diagonalized by per- 
forming three subsequent Bogolyubov transformations. 
First 



/?2 P t = 02 P a cos r? - (3 2pb sin rj 
fapl = fcpa sin rj + (3 2p b cos 77, 



(4) 



with cost? = , sin = and A 2 = \J (Af) 2 + (A|) 2 
Then 



ip1 



-op ' 



4 P sin ( 



/?2 P a = c a p sin 9 + d ap cos 9, 
fiipl = Cb P cos 9 - d bp sin 9, 
02 P b = c p sin -(- dbp cos 0, 



(5) 



with cos 6» = ^, sin6» = ^ an d A = -y/(A 2 ) z + (Af) 



And. finally. 



and 



aipt = p ap cos'i/' + e Q pSin?/>, 

"lpj. = Pbpcostp + eppSin-0, 

Cap = eapCOST/; - p ap sin -0, 

Cb p = e bp cos V> - Pbp sin ?/;, 



(6) 



(7) 



where if> = § + tp, cos 2 sin 2 ip 



1± 



V( £ p) 2 + A2 



and £ p = | [ep 1 — e£] The resulting quadratic Hamilto- 
nian becomes 



H eff _ /3 d f A + 

n 2 / t fc p"ttp""P^ 



XI S p ( e Ip e "P + ApVbp - e\ p e b p - pl p Pap ) , (8) 

p 



where E p = ±y (e p ) + |A| 2 . The self-consistent equa- 
tion for the gap reduces to 



1 



p a/( £p )^ + A 2 



(9) 



Two key observations follow from these results. First, 
self-consistency equation Eq.® sets the value of the to- 
tal order parameter A 2 + A 2 but does not specify what 
Ai and A 2 are. The implication is that, at this level, the 
OP manifold can be viewed as a 6— component vector 
(three components of Ai and 3 of A 2 ), SDW ordering is a 
spontaneous breaking of 0(6) symmetry, and an ordered 
state has five Goldstone modes. The degenerate ground 
state manifold is composed of two-sublattice states with 
antifcrromagnctic order along diagonals; (0,7r) or (jr, 0) 
states, for which Ai = or A 2 = 0, are just two of many 
possibilities (see Fig. [2|. Second, a linear combination of 
original electronic operators described by d a , P decouples 



(a) / X X X / / / 

X X X X X X X X 



/ X X X 

(c) \ X V X 

X \ X \ 

\ X \ X 

X \ X \ 



(d) 



X X X X 
X x x X 

X 



X 



X 



x^ 



X 



X 



x^ 



FIG. 2: (color online) Various SDW spin configurations de- 
scribed by A ie iQlR + A 2 e iQ2R . For the model of Eq. ©.only 
A? + Al is fixed. Panel (a) - Ai = 0, panel (b) - A 2 = 0, 
panel (c) - Ai _L A2, and panel (d) - Ai = A 2 



from the SDW procedure. Then, even in case of per- 
fect nesting, when sp" 1 = — £p, the system still remains 
a metal in the SDW phase - excitations described by 
e a ,p and p a , P operators become gapped, but excitations 
described by d a;P operators remain gapless. 

Note that the ground state degeneracy is even larger 
than in the J\ — J 2 model of localized spins - not only the 
angle between the two sublattices can be arbitrary but 
also the magnitudes of the ordered moments in the two 
sublattices can be different. To confirm this result, we 
computed the ground state energy E gr in the mean-field 
approximation and indeed found that it depends only 
on A 2 + A 2 . We also extended this analysis to include 
a — [3 interactions with momentum transfer (0,vr) and 
(it, 0) (£/ 2 terms in the terminology of Rebi2) and still 
found the same degeneracy. 

Note by passing that the momentum integration in the 
self-consistency condition, Eq. @ is not restricted to 
the FS, and a finite SDW order parameter A appears 
even when FS disappear. This is due to the fact that a 
particle-hole bubble made out of a and (3 fermions actu- 
ally behaves as a particle-particle bubble because of sign 
difference of fcrmionic dispersions of a and (3 fermions. 
The consequence is that the SDW ordered moment does 
not scale with the (small) sizes of the hole and electrons 
FSs. In other words, SDW order is due to dispersion 
nesting of a and (3 fermions, but not due to FS nesting. 



C. Selection of the SDW order 

1. Interaction between electron pockets 

We now add the interactions between the two elec- 
tron pockets and verify whether they lift the degener- 
acy. These interaction do not contribute to the quadratic 
form, but they do contribute to the A-dcpcndcnce of 



6 




the ground state energy. It has been argued^, based 
on the transformation of the underlying orbital model 
into a band model, that the interactions between the two 
electron pockets are not particularly small and must be 
included into the theory. 



FIG. 3: (color online) Energy dispersions in the SDW state 
for full nesting: (a) - one electron band and one hole band, 
(b) - one electron band and two hole bands. Dashed lines in 

(a) are the dispersions without SDW order. Thee are four possible p _ & interact i onS ; 



Hf -Ue^ 0ip 3 aPlpia'P^T>2<y'Plpia + U 7 ^ Ajpacr Plp 4 a' fcp2<7' PlpiO 

+ Y E [^aa/^/WApia + h.c] + Y E [E /^aa/^'/W/W + A^'/W/W] (10) 
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(we used the terminology consistent with Ref p^ 2 -) . It is Applying the sequence of Bogolyubov transformations 
natural to assume that all interactions are repulsive, i.e., and taking (•••), we obtain the contribution to the 
all Ui > 0. ground state energy from various terms in Eq. {TT 
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^E [E /3 lP3<J^lp4T'/ 3l P2^'/ 3 lPi^ +/ 3 2p 3 a^2p40-'/ 3 2p 2CT '/32p 1 c 



-AA' 



|A!| 2 |A 2 
A 4 



(11) 



r 



where dots stand for the terms which depend only on A 2 
and do not break a degeneracy, and A — (c — <£), where 



= ( c ! P c &p) = 

= (dLd b 



= M 



2N Sp I 1 



1 - 



and 
The 



l b P u b P I = (4p^t>p) = 2JV Sp 

quantity A is zero in the normal state, but has a finite 
value in the SDW state. Combining all contributions, we 
obtain 



TT<' J ' 



2A 2 [{U 6 + U 8 - 

(Ai • A 2 

AA 2 U 7 ± 



U 7 

2 



lAi 



A 4 



A 4 



(12) 



We 
A x • A 2 



see that E^. 

2 



depends on |Ai| 2 |A 2 | 2 and on 
i.e., it is sensitive to both, relative values 

and relative directions of Ai and A 2 . When all interac- 
tions are of equal strength, the first term vanishes, and 
the last term favors A\ _L A 2 . In this situation, the 0(6) 
degeneracy is broken, but only down to 0(3) x 0(2), i.e., 
the magnitude of the order parameter at each site is now 

the same because ^A\ + A 2 ^j = (&i — A 2 ^j , but the 
angle between the directions of the SDW order in the two 
sublattices (i.e., between Ai + A 2 and Ai — A 2 ) is still 
arbitrary. This is exactly the same situation as in the 
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classical J\ — J 2 model. However, once Us + Ug — Uj — U4 
is nonzero, the degeneracy is broken down to a conven- 
tional 0(3) already at the mean-field level. Because U4 
is reduced and even changes sign under RG — , while 
other Ui do not flow, the most likely situation is that 
U6 + U% — U7 — U4 > 0, in which case E e g % is minimized 

when either Ai = 0, or A2 = 0, i.e., SDW order is either 
(0,7r) or (n, 0). This is exactly the same SDW order as 
observed in the experiments. If Uq + U$ — U7 — U4 was 
negative, E g * would be minimized when |Ai| = IA2I, in 
which case the SDW OPs of the two sublattices would 
align orthogonal to each other. The spin configuration 
for such state is shown in Fig. [He). Such orthogonal 
spin configuration has been found in the analysis of spin 
ordering in the two-orbital model.— 



2. Deviations from, perfect nesting 

Consider next what happens when we also include into 
consideration the fact that electronic pockets are actually 
ellipses rather than circles, i.e., the effective masses m x 
and m y are not equal, and e^ 1 ^ e^ 2 . To continue with 
the analytical analysis, we assume that the ellipticity is 
small, introduce m x = (1+S)m and m y = (1— 8)m, where 
S << 1, and compute the correction to the ground state 
energy to second order in 8. Performing the same set of 
transformations as before, we find that, for a non-zero 5, 
Eq.© has to be supplemented by 

H P = 25 H P -^ L i cos2d t cos2 ^P e °p 

p 

+ sin 2 ipp\ v Pa V - di p d ap - sin -0 cos -0 (pi p e ap + h.c.)] 
+ sin 28 [cos-0 {p\pd ap + h.c.) - sin-0 (pt p d a p + h.c.)] } 

(13) 

where the angles 9 and ip arc defined in the same way 
as before, and the overall factor of 2 accounts for spin 
degeneracy. For a perfect nesting e Q [p a ] states are all 
empty [occupied], and for small 5 << A ellipticity does 
not change this. 

From Eq. (|13[) we obtain two contributions to E gr of 
order S 2 . One comes from virtual transitions to non- 
occupied states and is negative. Another comes from the 
change of the dispersion of the ungapped e^ 1 in the pres- 
ence of ellipticity, and is positive. The negative contribu- 
tion comes from non-diagonal terms in (|13p taken to sec- 
ond order. Applying a standard second order quantum- 
mechanical perturbation theory we obtain 




S 2 sin 2 26» x 
2 



A 



2 r 



2E r , 



(2Ep 



where e p and E p are defined after Eq. ([8]) , and dots stand 
for the terms which do not depend on 9 and do not break 
the degeneracy. We remind that cos 8 = |Ai|/A and 
sin6l = |A 2 |/A, so that sin 2 (26>) = 4|Ai| 2 |A 2 | 2 /A 4 . Re- 
placing the sum in Eq. (|14p by the integral and re-scaling, 
we obtain for the energy per unit area 



j^ja,ellipt 



-8 2 sin 2 28 F ( 4- 

16tt VA 



(15) 



where 



F(x) 



(y + x) 2 dy 

V1FT1 



(16) 



At large x, i.e., at small A/fi, expected within the itin- 
erant description, F(x -> 00) « 1, and Egf lhpt becomes 



E g ' r 



-S 2 sin 2 26> 



16tt 



(17) 



Another contribution to E gr comes from the diagonal 
term in Eq. (|13[) and is related to S— induced change in 
the dispersion of d— fermions, which arc not gapped by 
SDW. Adding -5(p 2 x -p 2 y )/m term from Eq.(Tr3]) to the 
dispersion of a d— fcrmion and evaluating the energy of 
the occupied states, we obtain after simple algebra that 
a linear term in S is canceled out, but S 2 term is finite 
and yields 



j^b,ellipt 



: 26» 



m/i 
~~8t7 



(18) 



Comparing the two contributions, we find that E^f 1 ^ 1 
is two times larger than Egf lhpt . Adding the two terms 
and expressing sin 2 28 in terms of Ai and A 2 , we obtain 



j^ellipt 



CjAiHAal 2 , C = S 



2 m M 
4ttA 4 



(19) 



(14) 



We see that C is positive, i.e., the correction due to el- 
lipticity of electron pockets breaks the degeneracy and 
selects either (0,7r) or (tt,0) state. This is again the 
same selection as one needs for consistency with the ex- 
periments. We found it quite remarkable that ellipticity 
introduces effective interaction between two SDW OPs 
which, for A << /i, leads to the same selection of the 
ground state SDW order as the direct interaction be- 
tween the two electron pockets. 

For larger A/fi, F(x) increases, eventually to F(x « 
1) = (l/2)(0.25/x) 2 , and E a g f Ulpt becomes larger than 
£b,eihpt _ rpj ic g jg n cnan g e f jrjeihpt 0CCurSj however, at 

quite large A > fi/2 which very likely is not realized in 
the pnictides. 
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3. Selection of the SDW order in the itinerant Ji — J2 
model. 



We next show that the selection of the (0, ir) and 
(%, 0) states in the itinerant model occurs only if fermion- 
fermion interactions are conventional charge-charge in- 
teractions rather than spin-spin interactions (the vertices 
contain spin S— functions rather than cr-matrices) . For 
briefty, we consider the itinerant J\ — J 2 model with spin- 
spin interaction and show that this model also possesses 
0(6) symmetry at the mean- field level, but the magnetic 
order selected by quartic terms is different from (0,7r). 

The itinerant J\ — J 2 model is described by 

H J *- J > =£s(p)S(-p)x 

[Jl (COS ^ + COSPy) + 2 J 2 COS p x COSp y )] (20) 

where 5(p) = (1/2) J2 P1 ^a^apOpi+pP^ where a are 
fermionic operators which can be either holes, when pi 
is close to (0,0), or electrons, when pi is close to (0, 7f) 
or (%, 0). The kinetic energy term is the same as in Eq. 

©• 

We focus on momenta p in (|2"0|) near (0,0), (0, tt), 
(tt, 0) and (tt, tt), re- write the Hamiltonian in terms of 
hole and electron operators, and introduce the same 
SDW vector OPs as before: Af = —Jsdw Sp( a ipT^ 1 pT) 



Substituting these coefficients into (|12|) . we obtain 



and A: 



:{x) 



-JsDwT,p( a ip^2 P ui)), where J S dw = 
(4J 2 — Ji). Keeping first only the interactions between 
a and /3 fcrmions and performing the same mean-field 
decoupling of the 4-fermion terms as in the previous 
section, we find that self-consistency equations for Ai 
and A 2 are again identical and only specify the value of 
A 2 = |Ax| 2 . 



|A 2 | 2 : 



Jsdw 
2N 



E 



1 



(21) 



The solution for A exists for J 2 > Ji/4, i.e. for large 
J 2 SDW order is antiferromagnetic along the diagonals. 
The condition J 2 > Ji/4 is similar to J 2 > Ji/2 for a 
classical J\ — J 2 model of localized spins (the conditions in 
itinerant and localized models do not have to be exactly 
the same, indeed). 

The degeneracy between different SDW states with the 
same A 2 is again broken once we include interactions be- 
tween electron pockets. Such interactions are generated 
by Ji and J 2 terms taken at momenta p = (0, 0) and 
(tt, tt). Rc-exprcssing S(p)S(~p) in terms of fcrmions 
we obtain the quartic interaction between (3\ and /3 2 
fermions in the same form as in Eq. Ijlip . with the coef- 
ficients 



U! 



Ji-3J 2 , £/ 7 s = 
3(Ji- J 2 ),UZ 



- Ji — 3J 2 , 
-3(Ji + J 2 ). 



E. 



J1-J2 



4Ji 



|Ai| 



IA2I 2 



A 4 



( Ji + 3 J 2 ) 



Ai • A 2 



A 4 



(23) 

where A > is defined after Eq. (fT2"|) . Comparing this 
form with Eq. (fT2|) we see that now the pre-factor for 

^Ai ■ A 2 ^j term is negative, i.e., the energy is lowered 



when Ai and A 2 are parallel. The terms |Ai| 2 |A 2 | 2 and 



Ai • A 2 ) arc then equal, and from ([23]) we obtain 



E J g r j2 



-12A 2 {J 2 - ,h) 



|Ai| 2 |A 2 | 2 
A 4 



(24) 



We see that the state with Ai = or A 2 = 0, which 
has (0, tt) or (tt, 0) order, is only favored in the range 
Ji/4 < J 2 < J 2 . For larger J 2 , the energy is minimized 
when Ai = ±A 2 . The corresponding SDW state has 
antiferromagnetic order for spins in one sublattice, but 
no SDW order for spins in the other sublattice (see Fig. 
[5£d)). Such state has been identified as one of possible 
candidates for the magnetic ground state in the generic 
analysis of the Landau theory for a two-component order 
parameter^ 

We see therefore that, when J 2 interaction dominates, 
the spin-spin interaction between the two electron bands 
selects different SDW order from the case when the in- 
teraction is in the charge channel. 



III. SDW ORDER IN FOUR-BAND MODEL 

So far we found that the stripe (tt, 0) or (0, tt) order is 
selected in the 3 band model (one hole and two electron 
FSs) with conventional charge interactions. We now add 
the second hole pocket and check how its inclusion affects 
the SDW order. 

As we said in the Introduction, the second hole FS 
is less coupled to electron FSs than the one that we al- 
ready included into the 3-band model. This is due to a 
combination of the two factors: the difference in the in- 
teractions U$dw an d the difference in the degree of the 
overlap with the elliptic electron FSs. We consider both 
factors. 

To begin, consider a model of two circular hole FSs and 
two circular electron FSs. Let's assume that all FSs are 
of the same size, but that there are two different SDW 



interactions between hole and electron bands - U. 

r{2} 



{1} 
SDW 



(22) 



for one hole band, nd Ug^ w for the other. We intro- 
duce four SDW OPs: An, Ai 2 , A 2 i, and A 22 , of which 
An and A 2 i are with momentum Qi, and Ai 2 and A 22 
are with momentum Q 2 . The OPs An and Ai 2 involve 
fermions from the first hole band, while A 2 i, and A 22 in- 
volve fcrmions from the second hole band. Without loss 
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of generality. An can be directed along z axis, and A\ 2 
in the xz plane, but the directions of A 21 and A 22 can 
be arbitrary in 3D space. To simplify the discussion, we 
assume that the SDW configuration is coplanar, and set 
An and A 2i to be along z-axis and Ai 2 and A 22 to be 
along x axis. In explicit form, we then have, by analogy 
with Eq. j3|), 

A U =A? = -^V5>I pT /3ipT> 
p 

p 

A 21 = Ai = -C/| 2 ZJ V^(4 p ^ lpT ) 
p 



A 22 = Af = -U™ w $>!p T /W- (25) 
p 



From now on the subindiccs 1 and 2 indicate SDW OPs 
associated with one or the other hole bands. 

As in previous section, we first consider only the inter- 
actions between hole and electron states which contribute 
to the SDW order (U\ and U 3 terms). Cecoupling 4- 
fcrmion terms using Eq. (|25[) . we obtain the quadratic 
Hamiltonian in the form H^ } f = H km + H^+^, 
where 



H 



kin 



E 



P, (7,2 — 1,2 



H i% = E HpT ( A ^ipT + A 2/3 2 pi) ~ a\ pl (Af(3 lpl - Af/3 2pT ) 



- E ["It ( A 3&pT + A 4/?2 P i) - 4 P i (AfApi ~ A^/3 2pT ) 



r{2} 



ft,.C. 

h.c. 



(26) 



(27) 



The part involving ct\ hole band and Af and A^ can be 
diagonalized in the same way as for the 3-band model, by 
introducing cos = Af/A^sin^ = Af/Ax where A x = 
\J (Af ) 2 + (Af) 2 and rotating j3i 2 into 



/3ipl = c b P cos 6» - d bp sin 6*, 

/?2 P | = c a p sin 6> + d ap cos 6 
/3i P T = c a P cos - d ap sin 6>, 
/?2 P t = -c bp sin( 



cos ( 



(28) 



Substituting this into (|26| we obtain 



#S = - Al E PipT c -p - <l c fcp] ( 29 ) 



and = H«™ c} + H«™ d} , where 



Tjkvn 



E/ £p [ a lpT ai PT a lpl ai Pi C ap c ap C op Cf, P 



(30) 



Tjkin 

E e P [ a 2 P T 



«2 P T + a 2pi a 2pi - 4 p ^ap - d f bp d bp 



(31) 



The part if/™ c i + involves hole ai operators and where i?i,p = ±\J (e p ) z + Af. This result is similar to 



electron c a ,6 operators: 



Tjkin i rr{2} 

E £ p [ a i P T ai PT + "lpi^lPl " C ip c ap - C^p^P 
P 



~ Al E [ a i P T Ca P ~ at ipl Cfc P 



(32) 



Eq. (|32p can be diagonalized by the same transformation 
as in Eqs. © and 0: 

"i P T = P ap cos ^ + e ap sin ip, 
Cap = e ap cos ip - pap sin ip, 
"ipl = Pb P cos ip + e bp sin ip, 
Cb P = ebp cos ip - p hp sin ^, (33) 



where, as before, ip 
1± 



+ ip and cos 2 ip, sin ?/> 



TjKtn | rr 



. Substituting, we obtain 

{2} _ 



E^ 



i.p 



e ap ea P + p\ p Pb P - e op e Dp - pl p p a p 



,(34) 



10 



Eq. ([8]). The self-consistent equation for the gap Ai is The remaining part of the Hamiltonian involves holes 

from a 2 band and reduces to 



U 



{i} 

SDW \ " 



2N ^ 



(35) 



P J £p +A( 



H 



9) + d a<p sin(0 - 



4 2pT \° a :P 



cos 



l 2pl V°a,p 



COS 



9) + 4, P sin(0 - 6>) 



(36) 



where we introduced cos 9 
A 2 



(Af) 2 . 



i|/A 2 ,sin( 



and 



This part of the Hamiltonian decouples from the one 
which we already diagonalizcd if we eliminate the terms 
with c— operators. This will be the case if we choose 
9 = tt/2 + 9. Then Eq. ([36]) becomes the quadratic form 
of hole «2 operators and electron d a ,b operators neither 
of which is present in the quadratic form of Eq. (f32|) . 
Diagonalizing the quadratic form in Eq. (|36[) . we obtain, 
in terms of new operators e a ,b and p a ,b. 



Tjkin 
• H {«2,d} 



H 



{2} 

a 2 /3 



E 



2,p 



clpGap 



PlpPbp 



PapPap 



,(37) 



where E 2 , p = ±y (e p ) 
tion for the gap A 2 is 

r{2} 



A|, The self-consistent equa- 



1 = 



U 



SDW 

2N 



E 



(38) 



We see that, with the choice oi 8 = tt/2 — 8, we decoupled 
the quadratic form into two parts, one for Ai, another 
for A 2 . One can easily check that for both parts it is en- 
ergetically favorable to have a state with non-zero A once 
self-consistent equations for Ai j2 have solutions. This im- 
plies that the energy is further reduced if, in addition to 
creating Af = Ai cos9 and Af = Ai sind, one also cre- 
ates A?; = — A 2 sin0 and Af = A 2 cos6'. Furthermore, 
for equal-size circular FSs, the solutions for Ai and A 2 
exist for any Ugp W and U^ w . The resulting fermionic 
excitations are fully gapped, i.e., within 4-band model 
of equal size circular hole and electron pockets, an SDW 
state is an insulator. 

We now verify whether such insulating state is consis- 
tent with the observed (0, tt) or (it, 0) order. Because the 
angle 6 disappears from the quadratic forms in Eq. (|32|) 
and Eq. (f3"6"|) (for 8 = ir/2 — 9), the ground state is again 
degenerate. The degenerate SDW OP manifold is given 
by 



5(R) oc n 2 (A lC os<9- A 2 sin(9)e 4QlR 



-Hx ( Ai sin 9 + A 2 cos ( 



D iQ 2 R 



(39) 



This degenerate set does contain (0, tt) and (tt, 0) states 
with only one ordering vector (either Qi or Q 2 ). These 
states are obtained if we set tan 8 = Ai/A 2 or tan 9 = 
— A2/A1. The issue then is whether such states are se- 
lected by other interactions in the same spirit as (0, w) or 
(tt, 0) orders were selected in the 3-band model. We argue 
that they are not. The argument is that (0, tt) and (tt, 0) 
ordered states in the 4-band model are obtained by choos- 
ing 9 such that Ai and A 2 have components with both 
momenta, Qi and Q 2 . The total spin component with 
either Qi or Q 2 vanishes because of cancelation between 
Ai and A 2 components. Indeed, for Ug 2 ^ w << u)^ w 9 
is close to zero or to tt/2, but it is not equal to either of 
these values. Meanwhile, when we considered lifting of 
the degeneracy in the 3-band model, we found that the 
interactions between /3i and /3 2 electrons and the ellip- 
ticity give rise to the 9 dependence of E gr in the form 
E gr (6) = E + Ex sin 2 (28) such that E gr (8) has min- 
ima at 8 exactly equal to or to tt/2. To verify what 

happens in the 4-band model, we extended those calcu- 

{2} 

lations to the case when Uq^ w is non-zero and analyzed 



-T2T 

the form of E gr (8) pcrturbativcly in Uc DW /U 



rfi> 



SDW I " SDW 



We 

skip the details of the calculations as they are similar to 
the ones for the 3 band model. We found that the min- 
ima in E gr (8) do not shift from 8 = and 8 = tt/2, at 



least at small U 



{2} 



w/^SD 



The consequence is that, 
for the 4 — band model, interactions which break the de- 
generacy of SDW manifold favor the states which are 
not (0, tt) or (tt,0) states. For example, when 8 = 0, 
S(R) ex Ain z e lQlR + A 2 n x e* Q2R . Such SDW state is 
a two-sublattice structure with equal magnitudes of the 
order parameters on the two sublattices, but with a non- 
collincar orientation of 

The outcome of this study is that (0, tt) or (tt, 0) order 
can only be preserved if A 2 is strictly zero, i.e., the second 
hole band is not involved in the SDW mixing. Only then 
interactions and ellipticity select (0, tt) or (tt, 0) states. 
Otherwise, the order necessary has both Qi and Q 2 com- 
ponents and the SDW OP has modulations along both x 
and y directions. 

This observation has an implication for the electronic 
structure. When A 2 = 0, one hole band and one electron 
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FIG. 5: (color online) The FS for the tight-binding model 
£ p = At' cos p x cos p y + 2t (cos p x + cos p y ) for t = (solid 
curves) and t = 0.5t' (blue curves), h and e refer to the hole 
and electron states occupying the corresponding parts of the 
BZ. Observe that the hole bands are now centered at (0, 0) 
and (n, n). 



FIG. 4: (color online) The changes to the magnetic structure 
introduced by the weak SDW coupling A2 between the second 
hole band and the electron bands. Panels (a) and (b) are for 
the FS topology of the pnictides, panels (c) and (d) are for 
the FS topology of t' - t model for t < t' . 



where e"k = e^ 1 cos 2 6 + e^ 2 sin 2 9 and £k 



cos 



e'l 1 sin 2 , 



The self-consistent equations for the two gaps 



band are not involved in the SDW mixing, i.e., the system 
remains a metal even when SDW mixing of the other two 
bands is strong. In other words, (0, 7r) and (7T, 0) SDW 
order necessary involves only one of two hole bands and 
only one of the two electron bands. The other hole and 
electron bands remain intact. 

As we said, for equal-size circular hole and electron 
pockets, there is a perfect nesting between both primary 
and secondary pairs of hole and electron states, and A 2 
is non-zero for any non-zero Ugjj w . This is, however, not 
the case when the two hole FSs are circles of unequal size, 
and the two electron FSs are ellipses. In this situation, re- 
doing the same calculations as above we obtain in general 
4 different SDW bands described by 



H S DW = J2J2 El ' 2 P\ e Up ea - b P + Pb,apPb,ap 

a,b p 



E 



3,4p 



( e a,hp e a 



bp + Pb,apPb,ap 



(40) 



Here 



Ei,2 P = - (et 1 + £k) ± 2 V & 1 - + 4|Ax| 2 (41) 



and 



E 



3,4p 



~ N 2 



±2V( £ k 2 ^) +4|A 2 | 2 , (42) 



1 = U 



{1} 

SDW 



1 = u 



{2} 
SDW 



E 

V 

E 



n (Ei p ) - n (E 2p ) 
K 1 -^k) 2 +4|A 1 p 
n(E 3p ) - n(E ip ) 

(^ 2 "^) 2 + 4|A 2 | 



(43) 
(44) 



Analyzing these equations, we found that, for a non- 
perfect nesting, SDW magnetism is a threshold phe- 
nomenon, i.e., to obtain a non-zero Ai and A 2 , the in- 
teractions Usdw must exceed the thresholds. We found 
that a hole band with a heavier mass, e.g., with a larger 
FS, is more strongly coupled to electron bands. This 
band then plays the role of the ct\ band in our the- 
ory. Once the corresponding Ug]^ w exceeds the threshold 

value Ucr \ the system develops an SDW order (Ai 7^ 0) 
This order, as we know, is a stripe (0, ir) or (ir, 0) order. 
For the smaller-size hole band (a 2 band in our termi- 
nology), SDW order with A 2 7^ emerges only when 



U 



{2} 
SDW 



exceeds a larger threshold Ucr > U ( 



Once 



this happens, SDW order acquires both Qi and Q 2 com- 
ponents. This is illustrated in FigfJJa)-(b). We see there- 
fore that the regime consistent with the experiments is 
Usdw > Ucr^ and Ug 2 ^ w < ■ This regime is quite 
realistic given that there is quite sizable difference be- 
tween the areas of the two hole pockets. 



12 



A. t-t' -U model 

In simple terms, our result that the stripe order is only 
consistent with A x ^ 0, A 2 is ultimately related to the 
geometry of the FS, namely to the fact that both hole FSs 
are located around the T point. Indeed, suppose that one 
hole and one electron FSs are mixed into an SDW state 
with the momenta between these FSs (say (Qi = (0,tt)). 
For strong interaction, the SDW excitations are gapped 
and effectively decouple from the other hole and electron 
bands. These two bands can also mix into an SDW state, 
however, because the second hole band is centered at 
(0,0), it is separated from the second electron band by 
Q2 = (tTjO), i.e., the second SDW order necessary has 
momentum (71", 0), and the stripe order gets broken. 

The SDW order would be different if the second FS was 
centered at (tt, it) because then the remaining hole and 
electron bands would be separated by the same Q as the 
first pair of bands, and the stripe configuration would not 
be broken by the emergence of the second SDW order. 
This is illustrated in Fig0Jc)-(d). 

Such FS geometry is realized in the half- filled t — t' — U 
model described by 

H = t c L c '+6i<? + *' X] C L C i+S2<r + U Y1 n ^ n n 

i.8\(T i,S2<7 i 

(45) 

where <5i and Si are distances to nearest and next-nearest 
neighbors, and n, = c\ a Ci a . In momentum space, the 
dispersion e p = 2t(cosp x + cosp y ) + At' cosp x cosp y — 11 
has maxima at (0, 0) and (n, tt) and isotropic quadratic 
hole-like dispersion around these points (with non-equal 
masses near (0,0) and (tt,tt), when t ^ 0), and minima 
at (0, tt) and (tt, 0) and elliptical electron-like dispersion 
near these points. We illustrate this in Fig. [5l Such band 
structure is topologically equivalent to the one shown in 
FigHc). 

At large U, the t — t' — U model reduces by standard 
manipulations to the J\ — J 2 Heisenberg spin model. The 
electronic states in this model are all gapped, and the 
SDW OP is degenerate at the mean-field level. We an- 
alyzed SDW order in the t — t' — U model for arbitrary 
U within the same computational scheme as before and 
found that (0, tt) and (tt, 0) states minimize the energy 
at the mean-field level because both Ai and A2 appear 
with the same momentum. 

It is instructive to consider the discrepancy between 
our 4-band model and t—t'—U model is some more detail. 
We remind that in our model, ellipticity of electron pock- 
ets and interactions between them select (0, tt) and (tt, 0) 
states already at the mean-field level. We show below 
that int—t' — U model there is no selection of a particular 
SDW order at the mean-field level even when t 7^ and 
electron dispersion is elliptical. In this situation, (0, tt) 
and (tt, 0) states remain degenerate with infinite number 
of other two-sublattice states. Beyond mean-field level, 
quantum fluctuations select (0, tt) or (tt, 0) order in the 
Ji — J 2 model. We haven't check the selection of the 



SDW order at smaller U, but it is likely that (0, tt) and 
(tt, 0) states are selected for all values of U. 

To demonstrate that mean-field SDW OP in t — t — 
U model remains degenerate, consider this model at 
small t and introduce two SDW order parameters Ai = 
-(U/2N)J2] s .( c L,a <T <xP c k+Qi,0) with momentum Qi, and 
A 2 = -(U/2N) J2k( c l, a (J af3Ck + Q 2 ,p) with momentum 
Q 2 . We introduce four fermionic operators with mo- 
menta near (0, 0), (0, tt), (tt, 0), and (it, tt), and re-express 
t — t' — U model as the model of holes and electrons. We 
skip computational details, which are not different from 
what we presented in previous Sections, and only cite 
the results. The extra terms in the quadratic form come 
from the interactions between holes and electrons. For 
t = 0, all 4 dispersions are equal (up to the overall sign), 
and the Hamiltonian is diagonalized in the same way as 
in Sec. IIIBI Not surprising, the ground state energy 
depends only on (Ai) 2 + (A 2 ) 2 , i.e., there exists a degen- 
erate manifold of SDW OPs. 

Consider next the effect of the interactions between the 
electron bands. There are four effective electron-electron 
interactions - the analogs of Uq, U7, Us and U4 in f| 10[) . 
They all originate from the same U term and have the 
same prefactors. According to Eq. (fT2|) . the extra term 
in the ground state energy from electron-electron inter- 
action only has (Ai • A 2 ) 2 part: 

(aVa,) 2 

E% = AA 2 U K A4 > (46) 

This term orders Ai and A 2 perpendicular to each other, 
what makes the magnitudes of the two sublattice OPs 
Ai + A 2 and Ai — A 2 equal. But in the absence of 
I Ai| 2 ■ |A 2 | 2 term, the angle between the two sublatticcs 
remains arbitrary. 

Consider next the effect of ellipticity of electron pock- 
ets. In our model, we remind, ellipticity gave rise to two 
contributions to E gr , both scale as |Ai| 2 • |A 2 | 2 . The 
two contributions were of opposite sign, but were not 
equal. In t — t' — U model, the situation is very sim- 
ilar - there are again two contributions to E gr : one is 
the second-order contribution from non-diagonal terms 
in the Hamiltonian, induced by t, another comes from 
the change of the dispersion of the diagonal terms. The 
contribution from non-diagonal terms is given by 

j?a,ellipt oj.2 |Ai| 2 |A 2 | 2 

V =8t a 3 

y (cos 2 Px + cos 2 p y ) A 2 + 
p [A 2 + 16(f) 2 cos 2 p x cos 2 p y f /2 

where, as before, the dots stand for the terms which de- 
pend only on A. The contribution from the change of 
the dispersion of the diagonal terms in the Hamiltonian 
is 
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E b,ellipt = 2 ^ 



E c l + £ c2 + £ fa + £ fb 



£ cl + £ fa 



£ c2 + £ fb 



(48) 



where 



e c i = At' cosp x cos p y + 2t (cosp x + cosp y ) , 
e c i = At cosp x cospy — 2t (cosp x + cosp y ) , 







//w\\ 


/ W\ 



Af-Al 



E/o = COS^ COSp a + 2t {cOSp x + COSp y ) I -^2 



£fb = —4t' cosp x cospy — 2t (cos p x + cospy) 



A? -A 2 



A 2 




(49) 



Substituting these energies into (|4"B")) and expanding in rj, 
we find 



-8t' 



|Ai| 2 |A 2 



A 4 



E 



(cos 2 p x 



COS' 



[A 2 



16(t') 2 COS 2 Pa; COS 2 Py 



3/2 



(50) 



rpa^ellipt 



Comparing (|4"7|) and ([SO)) , we find that 
j^MiKpt = o, i.e. , elfipticity of electron bands in t-t'-U 
model does not give rise to a selection of a particular 
SDW OP from the degenerate manifold. This result very 
likely holds for arbitrary t/t', as long as the ground state 
manifold consists of states with antiferromagnetic order 
along the diagonals. 



IV. ELECTRONIC STRUCTURE OF THE SDW 
STATE 

As we said earlier, when all four Fermi surfaces are 
circles of equal size, all 4 Fermi surfaces are mixed by 
the SDW and are gapped. However, as soon as the two 
hole pockets become non-equal, there is a range of the 
interactions when one hole pocket and one of the two 
electron pockets are mixed by SDW, but the other hole 
and electron pockets remain intact. As a consequence, 
in the folded BZ, the electronic structure contains one 
hole and one electron band mixed by SDW order and 
visible both near (0,0) and (vr, 7f), one hole band visi- 
ble only near (0,0), and one electron band visible only 
near (w, 7r). The only effect of SDW ordering on the non- 
mixed excitations is the A-dependence of the intensity 
via coherence factors. Besides, for moderately strong in- 
teractions, as in the pnictides, Ai is not too large, and 
the SDW-mixed dispersions form the new FSs^i. The 
net result is the presence of three bands near the T point, 



(e) 


(f) 


A-A 


/ w \ 



(0,0) 



FIG. 6: (color online) Calculated image of the electronic 
structure in the folded BZ in the SDW phase. Left panel 
- near F = (0,0), right panel - near M = In figures 

(a) and (b) we show the dispersion, in figures (c) and (d) the 
spectral function ImG(k, ui) = ^].lmGi, and in figures (e) 
and (f) ARPES intensity, i.e., the product IraG(k,u) xf(Ek). 
The direction of momenta in all figures is along k x = k y in 
the folded BZ. We directed SDW order parameter along z 
axis and set Af / 0, A 2 = A 3 = A 4 = ((0,7r) order). For 
definiteness we set Af = /i/4, \i\ — (3/2)/x, fi2 = (3/2)/x, 



and mi = m>2 ~ m y = 2m x 



1. For numerical purposes we 



added the small damping constant S = 0.04^. 



one of which (the hole band) is not modified by SDW, 
and another three bands close to the (ir, ir) point, one of 
which (the electron band) is not modified by SDW (see 
FigsU;a)-(b)). 

The spectral functions ImGi are peaked at quasi- 
particle energies but are also sensitive, via the SDW co- 
herence factors, to the interplay between the effective 
masses of the bare hole and electron dispersions. For ex- 
ample, for two equivalent hole pockets (toi = 7712) and 
m x < mi 



' u • 



the spectral weight of ARPES intensity 
^2 i riF(Ei ii ) x ImGi { u f{E) is a Fermi function), which 
we plot in FigsEle)-(f), is the largest for two hole-like 
bands centered around T-point, and for one electron band 
and hole "half-bands" (hole blades) around (ir, 7r)-point. 
This theoretical ARPES intensity is quite consistent with 
the experimental data from Refsi 28 ' 29 i 30 : the experiments 
also show two hole dispersions near T point, and one elec- 
tron band and two hole "half-bands" near (tt,tt). There 
is some evidence 2 ^ of the third hole dispersion near (0, 0), 
but this one likely emerges from the fifth band which we 
did not consider here. 
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E, 



(0,0) 



(Jl,7t) 



FIG. 7: (color online) Same as in Fig(B]but for a multi-domain 
sample with equal distribution of the domains with (tt, 0) and 
(0, tt) SDW orders. 



Our results are also consistent with recent observa- 
tions of "anisotropic Dirac cones" , 41 ' 42 Indeed, the bands 
mixed in the SDW state form small hole pockets near 
both T-and M-points of the BZ, and the dispersion 
around these pockets is almost linear close to the Fermi 
level and also anisotropic between k x and k y due to ini- 
tial anisotropy of the elliptic band involved in the SDW 
mixing. However, despite of a visual similarity with the 
Dirac cone, the dispersion is still quadratic near the top 
of the hole band. There is also an "accidental" Dirac 
cone in Fig. ([6]) near the T-point, at a momentum where 
the dispersion of the SDW mixed band intersects with 
that of the hole band which does not participate in the 
SDW mixing. This intersection need not to be at the FS 
but occurs close to the FS for the parameters chosen in 
Fig. ©. 

In Fig. [5] we assumed that a particular SDW order 
((0, 7r) in that figure) is the same in the whole sample. 
It is quite possible, though, that the system has domains 
with (tt, 0) and (0, tt) orders. This idea has been recently 
put forward in Refill. In Fig. [7] we show the disper- 
sions, the spectral functions, and the ARPES intensities 
for a multi-domain sample, assuming equal distributions 
of the domains with (7r,0) and (0,7r) orders. Comparing 
Figs. [5]and[71 we clearly sec that main effect of averaging 
over two different SDW orders fis the appearance of an 
dditional electron dispersion near the M-point. This sec- 
ond electron dispersion is consistent with recent ARPES 
measurements of the electronic structure in BaFe2As2 
by Shen and collaborators^ 1 - 



V. SUMMARY 

To summarize, in this paper we analyzed SDW order 
in the itinerant model for Fc-pnictidcs. We considered a 
model consisting of two hole bands centered at (0, 0) and 



two electron bands centered at (0, tt) and (tt, 0) points in 
the unfolded BZ (in the folded zone, this corresponds to 
two hole bands centered at (0, 0) and two electron bands 
centered at (7r,7r)). 

We assumed that one hole band is more strongly cou- 
pled to electron bands and considered first a 3 band 
model consisting of one hole and two electron bands. We 
found that SDW order in this model is highly degenerate 
for the case of a perfect nesting, and if we restrict with 
only the interactions between holes and electrons which 
give rise to a SDW instability. The degeneracy of the 
SDW order parameter is due to the fact that the energy 
gain to create SDW order parameters with momentum 
Qi = (0, tt) or Q 2 = (it, 0) is the same. The order param- 
eter manifold consists of states with a perfect antiferro- 
magnetic order along diagonals but with non-equal values 
and arbitrary directions of the OPs at nearest neighbors. 

We demonstrated that the degeneracy is broken once 
we include into consideration interactions between two 
electron pockets and cllipticity of these pockets. Both the 
interactions and the ellipticity select (0, tt) or (tt, 0) stripe 
phases with only one ordering momenta, either Qi or Q 2 . 
These stripe-ordered states are in agreement with the 
experiments on Fe As-based pnictides. The selection of 
either (0, tt) or (tt, 0) order leaves one electron FS intact, 
and the system remains a metal even when SDW order is 
strong, and the two bands involved in the SDW mixing 
are fully gapped. 

We argued that the selection of stripe states only oc- 
curs when the interaction between electron pockets is 
in the charge channel. We considered, as an example, 
the model with the interaction in the spin channel and 
showed that there electron-electron interactions select a 
different state. 

We next added the second hole band and argued that 
it mixes with the electron band which was left out of the 
SDW mixing in the 3 band model, once the corresponding 
coupling exceeds the critical value. We found, however, 
that this second SDW order necessary gives rise to a non- 
stripe spin configuration in which both ordering momenta 
Qi or Q 2 are present. The stripe (0,tt) or (tt,0) order is 
therefore preserved only if the interaction involving the 
second hole band is below the threshold, and the second 
SDW order does not develop. 

This in turn implies that, if the SDW order is (0,7r) 
or (tt, 0), as the experiments indicate, only one hole and 
one electron FSs are involved in the SDW mixing. Other 
two bands (one hole band and one electron band) are 
not involved and remain the same as in the normal state. 
As a result, the system remain a metal for any coupling 
strength. 

We argued that this peculiar requirement is related to 
the fact that both hole FSs are centered at (0,0). If 
one of two hole FSs was instead centered at (7r,7r), the 
second SDW order would have the same momentum as 
the first one (either Qi or Q 2 ), and the stripe order would 
survive. To illustrate this point, we considered a half- 
filled t — t U model in which the two hole bands are 



15 



centered at (0,0) and (%, 71"). In the large U limit, the 
model reduces to Ji — J 2 model of localized spins. We 
argued that, in this model, the degeneracy of the found 
state SDW order parameter manifold is not broken by 
either ellipticity or electron-electron interactions, and the 
degenerate manifold includes (0, w) and (ir, 0) states even 
when all four bands are involved in the SDW mixing. The 
degeneracy is broken by fluctuations beyond mean-field, 
which likely select the stripe order at any U. Because 
all 4 bands are involved, the stripe-ordered state in the 
t — t' — U modelit is a metal at small U and an insulator 
at large U . 

We analyzed the electronic structure for parameters 
relevant to the pnictides, for which SDW order is moder- 
ate, and the two FSs involved in the SDW mixing are only 
partly gapped. We found three bands near k = (0, 0) and 
three bands near k = (tt,tv) in the folded BZ, and more 
FS crossings than in the paramagnetic state. We calcu- 
lated ARPES intensity and found a number of features 
consistent with the data. In particular, we found "Dirac 
points" in the dispersion near (0, 0) and an electron band 
and two hole "half-bands" (hole "blades") near (tv,tt). 

We believe that the good agreement with ARPES ex- 
periments is a strong argument in favor of the itinerant 
scenario for the ferropnictides. We emphasize that the 
itinerant scenario does not imply that the system must 



be in a weak coupling regime. Interactions in the pnic- 
tides are moderately strong and, quite possibly, give rise 
to some redistribution of the spectral weight up to high 
energie s 9 ' 25 ? 43 . The only requirement for the applicability 
of our itinerant approach is the existence of a substantial 
spectral weight at low-energies, where the system behaves 
as an interacting Fermi liquid. 
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